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Abstract 

Extensive Air Showers are complex macroscopic objects initiated by single ultra-high energy particles. 
They are the result of millions of high energy reactions in the atmosphere and can be described as the 
superposition of hadronic and electromagnetic cascades. The hadronic cascade is the air shower backbone, 
and it is mainly made of pions. Decays of neutral pions initiate electromagnetic cascades, while the decays of 
charged pions produce muons which leave the hadronic core and travel many kilometers almost unaffected. 
Muons are smoking guns of the hadronic cascade: the energy, transverse momentum, spatial distribution 
and depth of production are key to reconstruct the history of the air shower. In this work, we overview the 
phenomenology of muons on the air shower and its relation to the hadronic cascade. We briefly review the 
experimental efforts to analyze muons within air showers and discuss possible paths to use this information. 



1 Introduction 

Our understanding of high energy physics is sup- 
ported by experiments up to the TeV scale. Beyond 
such high energy frontier we must rely on extrapo- 
lations of our theories in terrains which might hide 
unexpected phenomena. The only direct processes 
surpassing these high energies and which might rep- 
resent a challenge in particle physics are the reac- 
tions initiated by Ultra High Energy Cosmic Rays 
(UHECR) high up in the atmosphere. They reach up 
to ~ 10 20 eV in the lab system, which corresponds to 
a center of mass energy of about ~ 400 TeV. 

The origin and nature of UHECR remains a mys- 
tery. Our current understanding says that the vast 
majority of these particles are hadronic (atomic nu- 
clei [TJ [2]), excluding neutrinos [3] and photons [I]. 
The final solution to the UHECR puzzle must put to- 
gether different pieces: the astrophysical mechanisms 



that allow the acceleration to such gigantic energies, 
the propagation through the intergalactic space filled 
with magnetic fields, and last, the subject of this 
paper, the interaction with the Earth's atmosphere, 
which creates Extensive Air Showers (EAS). 

EAS encode the information of the primary among 
millions of secondaries by means of high energy inter- 
actions which lie on kinetic regions never accessed by 
experiments before. Among all secondaries, muons 
can travel many kilometers from the hadronic back- 
bone almost unaffected, carrying valuable informa- 
tion. Understanding this information is key to break 
the degeneracy between the uncertainties on the ex- 
trapolation of the hadronic interaction models to the 
highest energies and the composition of the UHECR 
beam. 

This paper is organized as follows: In section [2] we 
overview certain aspects of air showers and discuss 
the energy balance between the hadronic and elec- 



tromagnetic cascade. In section [3] we discuss how 
muons are produced in the hadronic cascade inher- 
iting valuable information from it. In section [4] we 
illustrate the ground distributions of muons. In sec- 
tion [5] we briefly discuss the experimental efforts to 
use this information. In section [6] we conclude. 

2 Extensive Air Shower dy- 
namics 

Extensive Air Showers are complex phenomenon ini- 
tiated by a single particle with an enormous energy. 
The collision with an air nucleus generates typically 
thousand of secondaries, which can interact again, 
creating a multiplicative process which is referred as 
cascade, and that can reach up to 10 11 particles at 
ground level for 10 20 eV showers. 

Depending on the kind of particles driving the mul- 
tiplicative process, there are two main subtypes of 
cascades. The ones initiated and driven by photons 
or electrons, and the ones originated and driven by 
hadrons. 

The study of the cascade can be done by means of 
the cascade equations, assuming some simplifications, 
or by means of full Monte Carlo simulations that in- 
clude many important details difficult to account for 
otherwise. On the other hand, Heitler models offer 
a simplified version of the main multiplicative pro- 
cess of a cascade and serves to qualitatively under- 
stand the most important features, giving approxi- 
mated values for relevant variables of the cascade. 
See for instance [5] for more details on the hadronic 
and EM cascade. 

2.1 The electromagnetic and the 
hadronic cascades 

When a high energy photon is injected into mat- 
ter, the most likely process to occur is an electron- 
positron pair production. Each of the new particles 
suffers bremsstrahlung, producing new photons. This 
multiplicative process repeats itself n times originat- 
ing the so called electromagnetic (EM) cascade. The 
total number of particles grows as 2™ . The energy of 



secondaries decreases as E = ^ to eventually reach 
the so called critical energy (E c ~ 80 MeV) at which 
electrons are more likely to lose their energy through 
ionization. At this point the cascade reaches the max- 
imum. After that, the multiplicative process stops, 
and the number of particles declines. The EM cas- 
cade practically keeps all the energy flowing within 
the EM channel, and does not leak into the hadronic 
cascade except for a small fraction by photopion pro- 
duction. 
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Figure 1: Energy fraction evolution with generation 
n. 

On a hadronic reaction at high energies, ~ 80% of 
the produced particles are pions, (7r + , ir~ and 7r°) in a 
~ 1:1:1 ratio, and ~ 8% are kaons, (K^, K s , K + and 
K~) also with a ~ 1:1:1:1 ratio. Neutrons/protons 
are produced with an overall probability of ~ 4-5%, 
and the rest is shared among other particles at the 
subpercent level as given by QGSJET-II.03 00. 

Neutral pions feed the EM cascade almost imme- 
diately, whereas charged pions either interact, sus- 
taining the hadronic cascade, or decay into muons 
(99.988%). In the same way, kaons interact feeding 
the hadronic cascade untill they reach their critical 
energy, which is of the same order of magnitude com- 
pared to pions. K s has a shorter lifetime (cr = 4 
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cm) compared to the rest of kaons (ct ~ few meters), 
which implies a higher probability of decay before in- 
teracting. A few hadronic generations after the first 
interaction, Kg decays, 31% of the times into 7r°7r°, 
and 69% into ir + n~. This means that kaons go from 
a ~ 0% contribution to the EM cascade in the first 
generations up to ~ 8% in higher generations, com- 
pared to the steady ~ 33% contribution of pions to 
the EM cascade. Finally, neutrons and protons keep 
interacting hadronically with no direct feeding into 
the EM cascade. 

The most relevant features of hadronic shower can 
also be approximately described by a Heitler model. 
After each hadronic generation n, there are cre- 
ated TO particles which subdivide in two main cate- 
gories: those which continue to feed the hadronic cas- 
cade, and those which feed the EM cascade, leaving 
the hadronic channel. They typically correspond to 
charged and neutral pions, in a 2 /3m and l/3m pro- 
portion. Thus, total number grows with the hadronic 
generation as (2/3m)™ whereas the energy decreases 
as E 71 = The energy fraction / carried by the 
sum of all charged pions in generation n to the total 
shower energy Eq is 

That is, in each generation, the energy carried by 
charged pions ^ IS reduced by a factor | . In a 
more realistic approach, we can include an effective 
factor k £ [0, 1] that modifies the amount of energy 
flowing to the EM cascade through n° decay as: 

/=(i-J«) n CO 

k can account for different aspects of the hadronic 
reactions. For instance, if a leading baryon takes 
(1 — k)Eq, k accounts for the inelasticity, being the 
fraction of energy going into pion production, and 
therefore \k,E goes into the EM channel, as ex- 
plained in [5], There might be other mechanisms 
that could effectively reduce the feeding to the EM 
channel, for instance, increasing the amount of kaon 
production [6j. 



2.2 The energy balance between cas- 
cades 

The energy share between both cascades evolves with 
the hadronic generation as showed in Fig. [T] In the 
beginning all the energy is in the hadronic sector. 
After 3 generations, (k — 1) 70% of the energy has 
been transfered to the EM sector. This means that 
the evolution of the EM cascade is rapidly decoupled 
from the hadronic cascade. Also shown is the case 
for k = 0.5, where the transfer from the hadronic 
to the electromagnetic cascade is slower. The energy 
balance affects the longitudinal developement and the 
muon content of the shower, see for instance [6]. 

The factor k can change with the energy of the 
hadronic reaction, and thus change with the gener- 
ation n. The energy at which the first and second 
generation reactions occur might be out of reach of 
the current man made accelerators. Fig. [I] also shows 
a case where k changes from a value K = 0.2to« = l 
after the first generation. It can be seen how the 
energy balance of the whole shower is affected. 

3 The production of muons in 
EAS 

Most muons in the shower come from the decay of 
pions, which are 10 times more numerous than kaons. 
Kaon decay can lead directly to muons (20%) or to 
charged pions (40%). 

Simple kinematics shows that the maximum trans- 
verse momentum p t that muons can obtain is just the 
center of mass momentum of the outgoing particles, 
which is 29.8 MeV. Given that the total momentum 
of the parent particles is of the order of a few tens 
of GeV, the direction of motion hardly varies, with 
deviation angle Af?^ ~ 0.01°. 

The experimental data of hadronic collisions avail- 
able up to a few hundreds of GeV per nucleon in the 
center of mass show a p t distribution that decreases 
exponentially 2 JjHdp t ^ ex P( — q) wnere Q changes 
slowly with the energy of the collision and the rapid- 
ity region. Q is of the order of tenths of GeV/c, that 
compared to the muon maximum transverse momen- 
tum available from the pion decay (~ 0.03 GeV) gives 
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Figure 2: Left panel: average, median and 90% quantiles of the y-distribution for different depths. Right 
panel: total number of muons produced per gem -2 , h(X), for 50 proton showers at 10 19 eV and 60 deg. 



a 10% correction. This makes the pt distribution of 
the outgoing muons very similar to that of their par- 
ents. This is a very important feature responsible for 
many of the observed characteristics of the hadronic 
and muonic showers. 

In [8] , it was argued that the transverse position of 
the production of muons, thus of the parent mesons 
decay, is confined to a relatively narrow cylinder: as 
the angle with respect to the shower axis goes as 
, the average traveled distance before the 
^tCTtt, where tv and m w are the 
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pion decay is I 
lifetime and mass of the charged pions. The perpen- 
dicular distance to the shower axis of the pion decay 
is r n = / sin a = T^2Q/m v ~ 22 m. 

Note that after each interaction n, the pt increases 
as pt ~ Q\fn. The outgoing angle goes as a geomet- 
rical progression with n as on ~ sina„ = (g^^~jn ■ 
The total outgoing angle a, is then dominated 

by the last interaction 



E* 



N 



{3/2N ch ) r 



(3) 



Fig. [2] (left panel), displays the y-coordinat^j] (the 
shower axis is at y=0) containing 50% and 90% of the 
production points as a function of the atmospheric 
depth. Also displayed is the average value, which is of 
tens of meters. This distance is small when compared 
to the distances involved in EAS experiments, which 
span from hundreds of meters to several kilometers 
in the perpendicular plane. For instance, the Pierre 
Auger Observatory has its tanks separated by 1.5 km 
[S]. Therefore, the position where the muon has been 
produced can be approximated by (0, 0, z), or simply 
z. 

Every cZA|^] along the shower axis, dN muons are 
produced within a given energy and transverse mo- 
mentum interval dEi and dp t . Their overall distribu- 
tion at production can be described in general with 
a 3-dimensional function, as: 



d 3 N 



dX dEi dept 



F{X,E,,cp t ) 



(4) 



x We use a system of coordinates (cylindrical (r, £, 2) or 
Cartesian (x,y,z), as convenient) with the z-axis aligned with 
the shower axis and (0,0,0) being at ground level. 

2 X is used as being equivalent to 2 corresponds to the 
column-air along z, X = p(z')dz' 
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The projection into the X (or z) axis becomes 

h(X) = J F{X,Et,cpt)dEidcpt (5) 

and it is the so called total/true Muon Production 
Depth (Distance) distribution, or MPD-distribution 
for short. It does not depend on the observational 
conditions since it does not contain any propagation 
effects of muons through the atmosphere. A detailed 
study of its shape is done in [TO]. Notice that this 
is different from the MPD-distributions of detected 
muons at a given position on ground gy|(r,c)) which 
includes the effects of propagation, as it will be ex- 
plained later. This distribution is sometimes referred 
to as apparent MPD-distribution. 

The total number of muons produced in a shower 

is 

Wo = y h(X)dX (6) 

It should be noted that this number is intrinsically 
different from the number of surviving muons, which 
is affected by the fluctuations of the depth of the first 
interaction, and thus change the distance traveled by 
muons to the ground. Some of the techniques used 
by experiments like Auger [TJJ use a fixed distance to 
the shower core, so they can also be affected by the 
lateral spread of the parent mesons. 

Eq. [4] can be factorized and expressed as the prod- 
uct 

F(X,Ei,cp t ) = h(X) f x (Ei,cpt) (7) 

where the function f x {Ei, cp t ) = ^/jrw^^ becomes 
the normalized Ei and cpt distribution at a given pro- 
duction depth X. In the approximations made in 
[8j [121 E] , fx did not depend on X and it was fac- 
torized in 2 independent distributions on Ej and cp t ■ 
This allowed analytical approximations of the distri- 
butions at ground. In |14j we have included these 
correlations, improving the accuracy of the energy, 
production depth, and time distributions at ground, 
and allowing for a proper description of the muon 
lateral distribution at ground. 

The function h(X) tracks the longitudinal devel- 
opment of the hadronic cascade and represents the 
production rate of muons per gem -2 . Its shape and 
features are extensively discussed in [10]. The depth 



at which h(X) reaches the maximum is denoted as 
^max- ^£ax correlates with the first interaction point 
Xi which corresponds to the first interaction of the 
primary in the atmosphere and the start of the cas- 
cading process [TU]. The most important source of 
fluctuations in air showers corresponds to the fluctu- 
ations of Xi, which causes an overall displacement 
of the whole cascade at first approximation. The 
amount X 1 = X — A^ ax defines the amount of tra- 
versed matter with respect to the shower maximum. 
The distributions can be expressed in terms of A', 
where the most important source of fluctuations has 
been eliminated, and only the remaining effects are 
present. 

In Fig. [2] (right panel) h(X) is shown for a sample 
of 50 showers. The fluctuations on the normaliza- 
tion and on A^ ax are clearly observed. In [T?| it was 
shown that both the energy and the transverse mo- 
mentum show similar features when referred to the 
same distance to the shower maximum, X' . 

In [51 [121 H3] the muon spectrum at production was 
approximated by a power law, E~ 2 ' 6 , following the 
high energy tails of the pion production on hadronic 
reactions. A more accurate description of this the 
spectrum was done in [14] : at low energies the single 
power law clearly does not work and, in addition, the 
energy spectrum evolves with X' by becoming softer, 
and stabilizing the shape after the shower maximum. 
In Fig. [3| left panel, the average energy spectrum 
of all muons at production is displayed for proton 
showers at 10 19 eV in different X' layers. 

The transverse momentum distributions are re- 
sponsible for most of the lateral displacement of 
muons with respect to the shower axis. In [SI H2 H2] , 
the pt distributions were approximated by an unique 
function, dN/dp t — pt / Q 2 exp(—p t / Q) , independent 
of the energy of the muon and its production depth, 
primary mass and zenith angle. In [14) . we uncover 
in detail all the dependencies. As the shower evolves, 
the pt spectrum becomes softer (Fig. [3j left panel 
shows the evolution as a function of X 1 ). Besides 
this dependence on X' , the p t distributions also de- 
pend on the energy of the muons, as discussed in 
[T4"] . The low energy muons display a smaller p t , and 
at high energies, the p t distribution prefers higher p t 
values. We have found that the different correlations 
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Figure 3: Normalized average energy (left panel) and average pt (right panel) distribution of all muons at 
production for proton showers at 10 19 eV and 60 deg zenith angle simulated with QGSJET-II.03 at different 
X' layers. 



of the pt with Ei and X must be included into the 
model in order to properly predict the muon lateral 
distribution at ground. 

In |14| it is shown that there are mild dependencies 
of both the energy and pt distributions on the energy 
and zenith angle of the primary. In addition, the 
photon initiated showers display quite different dis- 
tributions due to the different nature of the processes 
that lead to the muon production, through photopion 
production. Proton and iron showers, and different 
hadronic models also display mild differences among 
them. 

4 Propagation and ground dis- 
tributions 

In [H| it was shown that a few simple considerations 
are enough to account for most of the features ob- 
served in the muon distributions at ground. 

Firstly, muons exit the shower axis with an angle a 
determined by the energy and transverse momentum 
of the muon at production (sin a = ^f). The polar 



angle is distributed symmetrically over 2tt. Once the 
muon is produced, the trajectory is extrapolated in a 
straight line to the ground, and the arrival time due 
to geometric path is calculated. Once the main tra- 
jectory is defined, the energy loss, decay probability, 
multiple scattering and effects of the magnetic field 
are accounted for and the impact point on ground 
and arrival time delay are corrected. 

Table 0] summarizes different effects for 5 GeV and 
10 GeV muons produced at z = 10 km and arriving 
at a distance from the core r = 1000 m. 

The most important propagation effects that shape 
the ground distributions are, in this order: geometry, 
decay and energy loss. The magnetic effects become 
more important in showers with zenith angle above 60 
degrees. On the other hand, the multiple scattering 
effects are negligible at distances to the core above 
100 m. 

4.1 The energy distribution 

The energy at ground Ef was analyzed as a func- 
tion of the impact point on ground (r, £) . Typically, 
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Energy at production 


5.0 GeV 


10.0 GeV 


Energy at ground 


3.0 GeV 


7.8 GeV 


Probability of survival 


0.67 


0.84 


Geometric delay 


165 ns 


165 ns 


Kinematic delay 


12 ns 


2.3 ns 


Geomagnetic delay 


0.04 ns 


0.01 ns 


MS time delay 


1.5 ns 


0.8 ns 


Geomagnetic lateral deviation 


83 m 


17 m 


MS lateral smearing 


~ 60 m 


~ 35 m 



Table 1: Summary of the different effects after propa- 
gation for a muon produced at z=10 km and arriving 
at r=1000 m at 60 deg zenith angle, and geomagnetic 
field strengh perpendicular to the shower axis B± = 
10 /iT (MS stands for Multiple Scattering). 
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the muon energy is not directly measured by cosmic 
ray detectors since it would require carpeting exten- 
sive areas with particle detectors like those used in 
accelerator experiments. Nevertheless, the spectrum 
of muons has an impact on other quantities that are 
measured by current air shower detector arrays, like 
the muon lateral distribution at ground, the arrival 
angle, and the arrival time delay. 

Fig. [4] displays the normalized energy spectra of a 
60 deg shower, at different distances from the shower 
core. The energy of muons decreases as ~ 1/r and 
increases with the zenith angle [H [TH [T3] , being the 
details determined by the pt, z and Ei distributions. 
Low energy muons dominate at large distances from 
the core. 

4.2 Apparent production depth distri- 
bution 

The shape of the production depth distribution of 
the detected muons, the apparent MPD-distribution, 
changes with the observation position. The angu- 
lar position of the observation point respect to the 
production point z, selects particular (Ei,p t ) regions 
which can be more or less populated. In addition, 
the propagation effects, specially the decay, modu- 
late the apparent MPD-distribution depending on the 
energy spectrum of muons and also the path trav- 
eled from production to ground, I. Fig. [5] displays 



Figure 4: Normalized energy spectrum of muons ar- 
riving at ground for a 60 deg shower at different dis- 
tances from the core as given by CORSIKA compared 
to the prediction of the model. 



the apparent MPD-distributions for a 40 deg shower 
at different distances from the core, where the dis- 
tortions introduced in the dN / dX\^ r ^ distributions 
when compared to h{X) can be clearly observed. 

The dN/dX\r r n distribution is never directly ob- 
served, but reconstructed from the arrival time or 
the arrival angle at ground. The correct inference of 
the total/true MPD-distribution, h(X), requires the 
knowledge of the exact dependence of dN/dX\^Q 
with the observation point coordinates and detec- 
tion energy threshold. dN/dX\( r Q explores different 
kinematic regions at production when reconstructed 
at different distances from the core. For instance, the 
algorithm proposed in [TJ] and [TS] requires the con- 
version of each dN/dX\/ r ^) observed in each station 
to an universal distribution in order to sum up the 
contributions of all detectors in a single shower. 

4.3 Time distributions 

The total time delay is the sum of four different con- 
tributions t = t g + t e + tB + t-MS where t g is the 
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Figure 5: Comparison of several apparent MPD- 
distributions, dN/dX\r r Q, for a 40 deg shower at dif- 
ferent distances from the core. The total/true MPD- 
distribution (h(X)) is also plotted for comparison. 
Normalizations are arbitrary. 



geometric delay, i e is the kinematic delay, ts is the 
contribution produced by the geomagnetic field, and 
finally s includes the delay due to multiple scatter- 
ing. Fig. [6j left panel, displays the different contri- 
butions to the total delay for 60 degrees zenith angle. 
At large distances from the core, the geometric delay 
is the most important. At distances typically from 
a few hundred meters to 1 km, the kinematic delay 
has a large impact. As we increase the zenith angle, 
the geometric delay looses importance relatively to 
the other contributions. At 500 m from the core, the 
geometric delay represents ~60% of the total. Fig. 
[6j right panel, displays the overall time distributions 
at 1300 m from the shower core for a 60 deg shower. 
Filled histograms show the contributions of different 
muon energies at ground. High energy muons ar- 
rive earlier at ground. This is so because they are 
produced higher up in the atmosphere, and therefore 
have less geometric delay, but also because they have 
less kinematic delay. 

The muon arrival time distributions can be used 



to extract relevant information. Far from the core, 
the time distributions are to a very good extent a 
one to one map of the apparent MPD-distributions. 
They can be determined by converting each muon 
time into a production distance, being the kinematic 
time a second order correction. Since the energy 
of each muon is typically not known, it is approxi- 
mated by the mean value, taken from the energy spec- 
trum at each observation point as it was explained in 
[5J HH H3] . The energy would also determine the pa- 
rameters of the multiple scattering delay distribution, 
although its concrete value follows a random distri- 
bution. The geomagnetic delay can take only two 
possible values depending on the charge of the muon. 
In general this technique will require a stringent r 
cut for those regions where the geometric delay is a 
large fraction of the total delay, in order to avoid dis- 
tortions of the reconstructed dN / dX\( r _Q. A more 
promising method consists in fitting the time distri- 
butions at once leaving a set of shape parameters on 
h{X) free. Close to the core, the geometric delay 
is not dominant and the arrival time is mostly de- 
termined by the energy of each muon. This opens to 
possibility to measure, or at least constrain, the shape 
of the muon energy spectrum. A global fit would also 
allow to extract parameters from the p t distributions. 



4.4 Muon lateral distribution at 
ground 

The number of muons per surface area unit is 

p(r, C) — rdrcu; • ^ s ^ was snown m 03] J l° w en ~ 
ergy muons have a major impact on the fine details 
of the muon lateral distribution at ground. 

In vertical showers the number of muons per sur- 
face area does not depend much on £. As we increase 
the zenith angle, asymmetries appear because of the 
different propagation effects, mainly decay and geom- 
etry. The effects of the magnetic field become impor- 
tant above 60 degrees, and they completely dominate 
the distributions at very inclined showers, typically 
between 80 and 90 degrees [IB]. Fig. [7] displays the 
muon density as a function of r for 3 different polar 
angles £ on a 70 deg shower. 
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Figure 6: Left panel: Different contributions to the total average time delay for a 60 degrees shower. Right 
panel: Comparison between the model and CORSIKA of the normalized time distributions for a 60 degree 
shower r = 1300 m distances from the core. The color histograms show the contribution of different energies. 



The shape of the ground distributions is fully deter- 
mined by the distributions at production, h(X) and 
fx(Ei,p t ). A change in the overall muon content 
of the shower, Ao, produces a change in the muon 
density at ground, and therefore in the normaliza- 
tion of all distributions. The other main source of 
fluctuations comes from the depth of the first inter- 
action, which directly affects h(X) by changing its 
maximum, X^ ax . The position of IJJ^ directly in- 
fluences all distributions at ground since it changes 
the total distance traveled by muons to ground. 

4.5 Average energy and transverse 
momentum distributions 

One of the main applications of the present model is 
to be used in a global fit to extract information on 
the total number of muons in the shower Ao, and the 
total/true production depth distribution, h(X), and 
its maximum, X^ ay ,. In order to do so, a fx{Ei,pt) 
distribution must be assumed. 

The energy and transverse momentum distribu- 
tions display more universal features when they are 



expressed in terms of X' — X — X^ ax , once the effects 
of the fluctuations induced by the first interaction 
point are removed. The average energy and trans- 
verse momentum distributions do not change when 
changing the energy of the primaries, whereas they 
show mild differences between proton and iron pri- 
maries, and between hadronic interaction models. 

If we substitute fx(Ei,cpt) of a given shower by 
an average over showers of the same hadronic interac- 
tion model, primary, and zenith angle, (fx>(Ei, cpt)), 
leaving only h(X) from the original shower, the 
ground density displays differences of about ~ 2% 
at 1000 m compared to the prediction if we used 
fx(Ei,cpt), whereas the rest of the ground distri- 
butions remained unchanged. It is thus possible to 
use an universal energy and traverse momentum dis- 
tribution that depends only on X' , where the po- 
sition of A^ ax is naturally accounted for through 
X = X' + X^ ax . 

The systematics of any concrete application, in- 
cluding a global fit, are to be studied and accounted 
for in each particular method and/or experimental 
setup. The effects of the choice of hadronic inter- 
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action model on (fx>{Ei,cp t )) might introduce some 
systematics that should be also accounted for. On 
the contrary, those differences might be used to con- 
strain fx' itself when compared to data which is very 
promising. One could also think of a method to 
experimentally constrain the energy and transverse 
momentum spectrum based on simultaneous obser- 
vations of the ground distributions in different condi- 
tions. For instance, the ground muon distributions of 
inclined showers contain valuable information about 
the energy spectrum due to the spectrographic effect 
of the geomagnetic field. 

5 Experimental efforts 

In this section I will illustrate some of the experimen- 
tal efforts to reconstruct the muon distributions with 
a few selected examples. 

KASCADE has recently published [T7] the ap- 
parent MPD-distributions for showers between E~ 
[10 15 , 10 17 - 7 ] and zenith angle [0,18] deg at distances 
to the core [40,80] m. KASCADE uses a combina- 
tion of different detectors which can separate the 
components of the shower, being possible to indi- 
vidually tag single muons. It also has a muon tele- 
scope, able to track the trajectory of the muon back 
to the shower axis and thus determine the produc- 
tion height. The back-tracking technique can be used 
in combination with the time-to-X technique in the 
Time- Track Complementarity method |18j , which is 
able to separate high energy muons from low energy 
muons, opening new observables. As a drawback, 
these heavily instrumented observatories are hardly 
scalable to the high energy end of the spectrum, at 
energies around 10 19 eV. 

The Pierre Auger observatory has recently pub- 
lished the maximum of the apparent-MPD using the 
time-to-X technique. Although the water Cherenkov 
tanks were not specifically designed to distinguish 
muons from electrons and photons, a fiducial cut 
can remove those stations close to the core and keep 
the muon richness sufficiently high to reconstruct the 
MPD-distributions. 

KASCADE-Grande and Auger have also published 
in |19j and [20] the number of muons as a function of 
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Figure 7: Muon lateral distribution at ground for 3 
different polar angles £ for a proton shower at 10 19 
eV and 70 deg zenith angle. 



the energy, in two different energy ranges. Auger sees 
an excess in the number of muons when compared to 
simulations. It is still unclear whether the number of 
muons measured by the two experiments match due 
to the gap region around 10 18 eV. 

The shower-to-shower distribution of the number 
of muons contains also valuable information that it 
is not yet fully exploited. The RMS of the number of 
muons adds valuable information to help break the 
degeneracy between hadronic models and composi- 
tion. To achieve this goal, it is important to enhance 
the muon capabilities at the highest energies and gain 
precision in the muon reconstruction. 

MARTA (Muon Auger RPC Tank Array) is one 
of the efforts in this direction. It envisages to add 
highly robust and autonomous RPC detectors [H] 
to the Auger Cherenkov tanks to enhance the muon 
capabilities. 

6 Conclusions 

Particle reactions beyond the energy achieved by 
man made accelerators are continuously happening 
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in the hadronic backbone of EAS. Muons are true 
smoking guns of the hadronic cascade. In fact, we 
have shown that the relation between the distribu- 
tions directly inherited from hadrons and the distri- 
butions observed at ground are well understood and 
that ground distributions can be used to reconstruct 
and constrain the distributions at production. We 
must study muons to enhance the sensitivity to the 
hadronic phenomena in the cascade. This is the path 
to solve the problem of composition of the UHECR 
and might uncover new particle physics phenomena 
at high energies. The future cosmic ray experiments 
above LHC energies need precise muon dedicated de- 
tectors added to the EAS arrays. 
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